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Abstract 

A new type of supersymmetric transformations of the coupled-channel radial Schrodinger equa- 
tion is introduced, which do not conserve the vanishing behavior of solutions at the origin. Contrary 
to usual transformations, these "non-conservative" transformations allow, in the presence of thresh- 
olds, the construction of potentials with coupled scattering matrices from uncoupled potentials. As 
an example, an exactly-solvable potential matrix is obtained which provides a very simple model 
of Feshbach-resonance phenomenon. 
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During the past twenty years, Darboux transformations of the one-dimensional 
SchrMinger equation, also know, as supersymmetric quantum mechanics g have evolved 
into a powerful tool to solve the inverse scattering problem for a given partial wave |2J, |3( . In- 
deed, this formalism allows both the deductive construction from scattering data of a unique 
interaction potential without bound state jj| and the construction of all potentials, phase- 
equivalent to that potential but displaying different bound spectra [5|]. When the initial 
potential is zero, the transformed potentials not only fit experimental data with high preci- 
sion; they are also exactly solvable and have compact analytical expressions (|, a striking 
feature of the method which reveals its efficiency. 

This success seems however restricted to the single- channel case: for coupled channels, 
super symmetric transformations were not able up to now to provide a complete solution of 
the scattering inverse problem, despite the early generalization of their algebraic formalism to 
coupled equations Q. Whereas the construction of phase-equivalent potentials has proved 
to be possible [sj], a deductive construction of potentials by inversion of coupled- channel 
scattering matrices is still missing. In the case of coupled channels with different thresholds, 
such an inversion even seems impossible, as a matter of principle: applying a supersymmetric 
transformation to a diagonal (and hence uncoupled) potential matrix may result in a coupled 
ootential matrix |9| but the corresponding scattering matrix will always remain diagonal 
lo| . This makes it impossible to invert realistic data with coupling, which correspond to 
non-diagonal scattering matrices, by supersymmetric transformations of simple uncoupled 
potentials like the zero potential. 

The need for coupled-channel inversion techniques is however strongly growing these days 
since many problems in physics rely on coupled-channel systems. In nuclear physics for in- 
stance, theoretical calculations of hypernuclei require baryon-baryon interactions, which 
naturally display coupling between channels corresponding to different hyperon-nucleon and 
hyperon-hyperon states. In the context of experimental Bose-Einstein condensation, atom- 
atom interactions are monitored with the help of magnetic Feshbach resonances Q], a phe- 
nomenon based on the coupling between different spin states displaying distinct thresholds 
in the presence of an external magnetic field 

The purpose of this Letter is to show that the limitations on coupled-channel supersym- 
metric transformations mentioned in Ref. do not hold for a category of transformations 
not used up to now. These transformations do not conserve boundary behaviors of solutions 
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of the Schrodinger equation at the origin, which is why we call them "non-conservative" . We 
establish general properties of these new transformations, compare them with usual trans- 
formations and finally apply them in the simplest possible situation: a single transformation 
of an initial vanishing potential. This leads to a Feshbach-resonance exactly-solvable model. 



Using definitions from coupled-channel scattering theory 
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3, 



we consider a set of iV 



coupled radial Schrodinger equations, that read in reduced units 

- ip"(k, r) + V(r)ip(k, r) = k 2 ip(k, r), (1) 

where prime means derivation with respect to the distance r between the two bodies. Poten- 
tial V is an N x N real symmetric matrix, supposed to be short ranged and bounded. The 
complex wave-number diagonal matrix k is defined by its elements ki = \JE — Aj, where E 
is the center-of-mass energy and = Ai<A2<...<An are the threshold energies of the 
iV channels. The solution ip is an N x iV solution matrix made of iV solution vectors. 

Equation being of second order, its general solution can be constructed by linear 
combination, with matrix multiplication on the right-hand side, of two solution matrices 
made of 2N linearly-independent vectors. In the following, we use two such pairs of solutions: 
Jost solution matrices f(±k,r), that behave asymptotically as pure exponentials, 

f(±k,r) ~ exp(±zA;r), (2) 

and a pair of solutions made of the regular solution matrix <p(k,r), which vanishes at the 
origin, and solution matrix rj(k,r), the derivative of which vanishes at the origin. We 
normalize these as (p'(k, 0) = rj(k, 0) = /, where I is the identity matrix. In terms of Jost 
solutions, the regular solution reads [l^ 

<p(k, r) = i [f(k, r)fc- 1 F(-A;) - f(—k, r)fc" 1 F(A;)] , (3) 

which defines the Jost matrix F(k) as 

F(k) = f(k,0), (4) 

where T means transposition. This is obtained by calculating, both at the origin and at 
infinity, the Wronskian W[ip(k,r), f(k,r)] = if T (k,r)f'(k,r) — (p lT (k,r)f(k,r), the value of 
which is independent of r, taking into account W[f(—k,r),f(k,r)] = 2ik. The Jost matrix 
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is a key quantity in scattering theory: bound (resp. resonant) states correspond to zeros of 
its determinant in the upper (resp. lower) k{ planes and the scattering matrix S reads 

S(k) = k- 1/2 F(-k)F~\k)k 1 ' 2 . (5) 

For that reason, the following study mostly concentrates on Jost-matrix properties. 

To perform a supersymmetric transformation of Eq. we follow Ref. Q and first 
factorize it as 



A + A- - k 



ip(k,r) = k 2 ip(k,r), (6) 



where = ±d/dr + U(r) are mutually-adjoint first-order differential operators, defined 
in terms of superpotential real symmetric matrix U. Equation (JBJ) is equivalent to Eq. (JTJ 
when 

U(r) = a'(r)a- l {r), (7) 

where a(r) is a solution matrix of Eq. at negative energy S, referred to as the factorization 
solution. This solution matrix is assumed to be real and invertible for all r; moreover, its self- 
Wronskian W(a, a) vanishes in order for U to be symmetric. Another important property of 
Eq. (0) is that multiplying the factorization solution on the right by an arbitrary constant 
regular matrix, which is equivalent to changing the factorization solution, does not affect 
the superpotential. This strongly reduces the number of possible transformations. Finally, 



the wave-number diagonal matrix k is defined by its positive elements K{ = \/A» — 8. 

We then apply A~ on the left to Eq. (jHJ), which leads to a new Schrodinger equation with 
potential matrix 

V(r) = V{r) - 2U'(r) = -V(r) - 2k 2 + 2U 2 {r) (8) 

and solutions 

4>(k,r) = A~ip(k,r). (9) 

The asymptotic behavior of this relation when if} is a Jost solution of the initial equation 
leads to the following expression for the Jost solution of the new equation 

f(k,r)=A-f(k,r)[U( 00 )-ik]- 1 . (10) 



As will be proved elsewhere, when all thresholds are distinct, U(oo) is a diagonal matrix 
with elements ±«i, the signs depending on the asymptotic behavior of factorization solution 
a. Let us stress that Eq. (jlUJl shows that supersymmetric transformations always transform 
a Jost solution into a Jost solution, up to a multiplicative factor; in particular, when / 
vanishes at infinity, / vanishes too. 

The situation is very different for the behavior at the origin: when applied to a regular 
solution of the initial equation, Eq. Q does not systematically lead to a regular solution 
of the new equation, which leads to the distinction between two types of transformations. 
We call "conservative" transformations those that transform a regular solution of the initial 
equation into a regular solution of the new equation: (p{k,r) oc A~ip(k, r). Most transfor- 
mations used up to now in the literature j?l Q belong to this category. Using Eqs. (J3J) and 
(111, these transformations can be shown to lead to simple modifications of the Jost matrix, 

nn 

similar to the single-channel case [16|] , 



In this equation, the upper (resp. lower) signs correspond to a diverging (resp. vanishing) 
at the origin, two cases which will be considered in detail elsewhere. For vanishing cr(0), 
Eq. m , combed with Eq . B, leads for instance to the scattering-matrix modification by 
supersymmetric transformations obtained in Ref. p|. 

For one-channel elastic scattering, Eq. (jlljl is most useful in the context of inverse scatter- 
ing problem: such transformations can easily be iterated, which leads to a Pade approximant 
of the Jost matrix of arbitrary order, used to fit experimental scattering data [3(. In the 
coupled-channel case with distinct thresholds, however, this Jost-matrix modification is very 
restrictive: since both k and U(oo) are diagonal matrices, an initial decoupled potential al- 
ways transforms into a potential with a decoupled Jost matrix, and hence a decoupled 
scattering matrix, as found in Ref. Q]. Consequently, these transformations are not able 
to fit data corresponding to non diagonal scattering matrices. 

In the present work, we consider a new category of supersymmetric transformations: those 
for which the factorization solution is a finite regular matrix at the origin. In this case, the 
value of the superpotential at the origin, £7(0), is also finite and can be fixed arbitrarily. 
This can be shown by expressing a as a linear combination of (p and rj satisfying Eq. (J2J): 




(11) 




(12) 
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where the assumption cr(0) = / does not reduce the generality of the superpotential. Apply- 
ing Eq. Q to a regular solution of the initial equation leads in this case to a solution of the 
new equation which does not vanish at the origin. This type of supersymmetric transfor- 
mations hence break boundary conditions and we propose to call them "non-conservative" . 
Similar transformations have already been used in Ref. Q| but for the very particular case 
of zero factorization energy in the single- channel case. Using Eqs. (j3J) and (|10p. one gets the 
modification of the Jost matrix for non- conservative transformations 

F(k) = [U{oo) - tk}- 1 [F(k)U{0) - f' T {k, 0)] . (13) 

This relation is clearly more complicated than for conservative transformations, with the 
drawback that its iteration might be difficult. However, it presents the enormous advantage 
that even when the initial potential is uncoupled, which means its Jost solution and matrix 
are diagonal, the transformed Jost matrix becomes coupled by choosing a non-diagonal U (0). 
This new type of supersymmetric transformation thus seems to solve the main drawback 
of supersymmetric quantum mechanics for coupled-channel systems with distinct thresholds 
stressed in Ref. 

Let us now illustrate these findings by applying the above formalism to an initial vanishing 
potential V = 0, for which F(k) = S(k) = I. The most general factorization solution matrix 
for a non- conservative transformation can be written in two equivalent ways, 

cr(r) = cosh(/tr) + sinh(fi;r)K _1 [/(0) 

= exp(nr)C + exp(— kt)D. (14) 

For N channels, the transformed potential, as obtained from Eqs. (0IHI), depends on 
Ki,...,Kjf (or equivalently on iV — 1 thresholds and one factorization energy) and on 
N(N + l)/2 arbitrary parameters appearing in U(0). It is exactly solvable and is equivalent 



to the potential derived by Cox |l8l . Il9j in his q = 1 case. However, the present derivation 
not only leads to a much simpler analytical expression for the potential [compare Eq. (4.7) 
of Ref. [3| to our Eqs. (JJJ) and (JHJ)] but also subsumes it: the restriction det A ^ of Ref. 
which is equivalent to detC ^ in our Eq. (JHJ), does not apply here. 
In the case det (7^0, one has U(oo) = k, as shown by Eqs. (J2J) and (|14jl . According to 
Eq. ()13j) . the Jost matrix of the new potential then reads 

F(k) — (k — ik)" 1 [U (0) - ik] . (15) 



For N = 2 channels for instance, the potential depends on five real parameters 

/ 




17(0) 




(16) 



and Eq. (JT5|l is equivalent to Eq. (5.1) of Ref. [laj . 

As stressed above, our formalism is also valid for det C 
U(oo) = diag(— Ki, K2), which leads to the Jost matrix 



0. In this case, one has 



/ 



F(k) 



i/3 



10 fc 2 +»«2 



(17) 



This new result cannot be obtained from Cox' formula. Moreover, for det C = det D = 0, 
which is equivalent to the particular choice of parameters 



"1,2 = ±\J (K1K2 - P 2 )(ki/k 2 ) ±1 («i«2 > (3 2 ), 



(18) 



further simplifications occur. One gets superpotential 



U{r) 



coshy 



-K\ sinh y i/KiKq 



(19) 



with y = (« 2 — Ki)r — arccoshy kik 2 /(3 2 and potential 



V 



2{k 2 - wi) 
cosh 2 y 



K\ sinh y 

i y/KiK2 sinh?/ — 



(20) 



where diagonal term Vn (resp. V22) is positive (resp. negative). The determinant of Jost 
matrix (|T7j) has two zeros: one in k 1R = K2/K1P — k 2 R = —J K1/K2P — 10:2, the 
other one in — /c* R , — A^r- Since these zeros lie in the lower-half fcj and in the upper-half &2 
complex planes, they correspond to a resonance in channel 1, only visible below threshold 



A 



The three remaining parameters Ki, k 2 and (3 can be expressed in terms of 



this resonance energy E^ and width T, as defined by k\ K = £r — iT/2. They read 



2k 12 
4/3 4 



El + P/4 + J(E R - A) 2 + T 2 /4 T A, 



£ R + ^El + T 2 /4 
x ( E n - A + J(E R - A) 2 + T 2 /4 



(21) 
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Let us stress that, in general, it is not trivial to calculate the zeros of the determinant of 
the Jost matrices (fTSJl and (fT7j). i.e., to calculate the location of their bound states and/or 
resonances. In the spirit of the inverse problem, it is however very useful to construct po- 
tential models with bound and resonant states corresponding to physical values. The above 
potential provides such a model: it is a simple, exactly-solvable, two-channel interaction 
with one controllable resonance. When En < A, this resonance is a Feshbach resonance, for 
which the model constitutes a good pedagogical example: when T = 0, all elements of V 
vanish, except for V22 which has then a bound state at energy E-&. This corresponds to a 



zero of the Jost-matrix determinant in km = v^r, &2R = t\/ A — Er. When T > 0, coupling 
turns this bound state into a resonance by moving the zero off the axes, counter-clockwise for 
k-i and clockwise for k\. For instance, for A = 10, £7r = 7 and T = 1, one gets the potential 
plotted in Fig. Q The corresponding scattering matrix, as defined by its eigenphases and 



mixing parameter 



□ 



19j, is represented in Fig. El where the Feshbach resonance can be 
seen on b\. Crossing threshold, which corresponds to going from positive imaginary k 2 to 
positive real fc 2 , produces an interesting cusp effect in 8\ for which the present model 
also constitutes an analytical example. The case E-& > A, though allowed, is less interesting 
from the physical point of view: when T = 0, coupling does not disappear and the zero of 



the Jost-matrix determinant lies in km = y/E-^, /c 2 r = — \/-Er — A, with no strong impact 
on the physical region ki > \/A, k 2 > 0. 

In conclusion, we have introduced a new type of supersymmetric transformations which 
are able to transform an uncoupled potential into a potential with a non-diagonal scattering 
matrix for non- vanishing thresholds (see the non-trivial behavior of e in Fig. EJ). This is 
a promising first step for coupled-channel inversion in the presence of thresholds. The 
simplest possible application of this formalism leads to an exactly-solvable potential which, 
for particular choices of parameters, provides a textbook example of the Feshbach-resonance 
phenomenon, with compact analytical expression both for the potential and its Jost matrix. 
Future research will try to extend this result to iteration of transformations, with the hope 
to get a supersymmetric-quantum-mechanics inversion technique as efficient in the coupled- 
channel case as it is in the single- channel case Q|. 

This text presents research results of the Belgian program P5/07 on interuniversity at- 
traction poles of the Belgian Federal Science Policy Office. B.F.S. is partially supported by 
RFBR grant 06-02-16719 and thanks the National Fund for Scientific Research, Belgium, 
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FIG. 1: (Color online) Two-channel potential matrix (|20|) as denned by parameters (|21|) for 
threshold energy A = 10, Feshbach-resonance energy £"r = 7 and width T = 1. 
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FIG. 2: (Color online) Eigenphases S\, 62 and mixing parameter e of the scattering matrix corre- 
sponding to the potential of Fig. UJ as defined by Eqs. ©, (fT7|l. ifTHl) and (f2T|) . 

for support during his stay in Brussels. 
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